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We study the propagation and localization of classical waves in one-dimensional disordered struc- 
tures composed of alternating layers of left- and right-handed materials (mixed stacks) and compare 
them with structures composed of different layers of the same material (homogeneous stacks). For 
weakly scattering layers, we have developed an effective analytical approach and have calculated the 
transmission length within a wide range of the input parameters. This enables us to describe, in a 
unified way, the localized and ballistic regimes as well as the crossover between them. When both 
refractive index and layer thickness of a mixed stack are random, the transmission length in the 
long-wave range of the localized regime exhibits a quadratic power wavelength dependence with dif- 
ferent coefficients of proportionality for mixed and homogeneous stacks. Moreover, the transmission 
length of a mixed stack differs from the reciprocal of the Lyapunov exponent of the corresponding 
infinite stack. In both the ballistic regime of a mixed stack and in the near long-wave region of a 
homogeneous stack, the transmission length of a realization is a strongly fluctuating quantity. In 
the far long-wave part of the ballistic region, the homogeneous stack becomes effectively uniform 
and the transmission length fluctuations are weaker. The crossover region from the localization to 
the ballistic regime is relatively narrow for both mixed and homogeneous stacks. In mixed stacks 
with only refractive-index disorder, Anderson localization at long wavelengths is substantially sup- 
pressed, with the localization length growing with wavelength much faster than for homogeneous 
stacks. The crossover region becomes essentially wider and transmission resonances appear only in 
much longer stacks. The effects of absorption on one-dimensional transport and localization have 
also been studied, both analytically and numerically. Speciflcally, it is shown that the crossover re- 
gion is particularly sensitive to losses, so that even small absorption noticeably suppresses frequency 
dependent oscillations in the transmission length. All theoretical predictions are in an excellent 
agreement with the results of numerical simulations. 



PACS numbers; 42.25.Dd,42.25.Fx 
I. INTRODUCTION 

Metamaterials are artificial structures having negative 
refractive indices for some wavelengths^. While natural 
materials having such properties are not known, it was 
the initial paper— that sought to realize artificial meta- 
materials which triggered the rapidly increasing interest 
in this topic. Over the past decade, the physical proper- 
ties of these structures, and their possible applications in 
modern optics and microelectronics, have received con- 
siderable attention (see e.g. Refs^"— ). The reasons for 
such interest are their unique physical properties, their 
ability to overcome the diffraction limit^'^, and their 
potential role in cloaking^, the suppression of sponta- 
neous emission rate^, and the enhancement of quantum 
interference^, etc. 

Until recently, most studies considered only ideal sys- 
tems and did not address the possible effects of disor- 
der. However, real metamaterials are always disordered, 
at least, in part, due to fabrication errors. Accordingly, 
the study of disordered metamaterials is not just an aca- 
demic question but is also relevant to their application. 
The first step in this direction was made in Ref:-*^ where 



it was shown that the presence of a single defect led to 
the appearance of a localized mode. A metamaterial with 
many point-like defects was considered in Ref.^^ where it 
was demonstrated that even weak microscopic disorder 
might lead to a substantial suppression of wave propaga- 
tion through a metamaterial over a wide frequency range. 

The next steps in this direction focused on the study 
of localization in metamaterials. Anderson localization 
is one of the most fundamental and fascinating phe- 
nomenon of the physics of disorder. Predicted in the 
seminal paper^^ for spin excitations, it was extended to 
the case of electrons and other one-particle excitations 
in sohds, as well as to electromagnetic waves (see, e.g., 
Refsiii^"— ), becoming a paradigm of modern physics. 

Anderson localization results from the interference of 
multiply scattered waves, manifesting itself in a most pro- 
nounced way in one-dimensional systemsliii^, in which 
all states become localized^^ so that the envelope of 
each state decays exponentially away from a randomly 
located localization center—. The rate of this decay is 
non-random and is called the Lyapunov exponent, 7, the 
reciprocal of which determines the size of the area of lo- 
calization. 
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In a finite, but sufficiently long, disordered sample, the 
localization manifests itself in the fact that the frequency 
dependent transmission amplitude is (typically) an ex- 
ponentially decreasing function of the sample size. The 
average of this decrement is a size-dependent quantity, 
whose inverse (i.e., reciprocal) is termed the transmis- 
sion length, 1^. In the limit as the sample becomes of 
infinite length, the decrement tends to a constant non- 
random value. The reciprocal of this value determines an- 
other characteristic spatial scale of the localized regime, 
which is the localization length, I. It is commonly ac- 
cepted in both the solid state physics and optical commu- 
nities, that the inverse of the Lyapunov exponent, 7~^, 
and the localization length, I are always equal. While 
this is true for media with a continuous spatial distri- 
bution of the random dielectric constant, in the case of 
randomly layered samples, the situation, as we show in 
this paper, is more complicated. In particular, the in- 
verse of the Lyapunov exponent by itself, calculated, for 
example, in Refi^, does not provide comprehensive in- 
formation about the transport properties of disordered 
media. Furthermore, it is unlikely that it can be mea- 
sured directly, at least in the optical regime. 

The first study of localization in metamaterials was 
presented in Refi^ where wave transmission through an 
alternating sequence of air layers and metamaterial lay- 
ers of random thicknesses was studied. Localized modes 
within the gap were observed and delocalized modes were 
revealed despite the one-dimensional nature of the model. 
A more general model of alternating sequences of right- 
and left-handed layers with random parameters was stud- 
ied in Refi^. There, it has been shown that in mixed 
stacks (M-stacks) with fluctuating refractive indices, lo- 
calization of low-frequency radiation was dramatically 
suppressed so that the localization length exceeded that 
for homogeneous stacks (H-stacks), composed solely of 
right- or left-handed slabs, by many orders of magnitude 
and scaling as Z oc or even higher powers of wave- 
length (in what follows we refer to this result as the 
A® anomaly), in contrast to the well-known dependence 
I (X A^ observed in H-stacks^. As noted in KeiM, a 
possible physical explanation of this is the suppression of 
phase accumulation in M-stacks, related to the opposite 
signs of the phase and group velocities in left- and right- 
handed layers. Scaling laws of the transmission through 
a similar mixed multilayered structure were studied in 
Refi^. There, it was shown that the spectrally averaged 
transmission in a frequency range around the fully trans- 
parent resonant mode decayed with the number of layers 
much more rapidly than in a homogeneous random slab. 
Localization in a disordered multilayered structure com- 
prising alternating random layers of two different left- 
handed materials was considered in Ref.^^, where it was 
shown that within the propagation gap, the localization 
length was shorter than the decay length in the underly- 
ing periodic structure, and the opposite of that observed 
in the corresponding random structure of right-handed 
layers. 



In this paper, we study the wave transmission through 
disordered M- and H-stacks of a finite size composed of 
a weakly scattering right- and left-handed layers. In the 
framework of the weak scattering approximation (WSA), 
we have developed a unified theoretical description of the 
transmission and localization lengths over a wide wave- 
length range, allowing us to explain the pronounced dif- 
ference in the transmission properties of M- and H-stacks 
at long wavelengths. 

When both refractive index and layer thickness of the 
mixed stack are random, the transmission length in the 
long wavelength part of the localized regime exhibits a 
quadratic power law dependence on wavelength with dif- 
ferent constants of proportionality for mixed and homo- 
geneous stacks. Moreover, in the localized regime, the 
transmission length of a mixed stack differs from the re- 
ciprocal of the Lyapunov exponent of the corresponding 
infinite stack (to the best of our knowledge, in all one- 
dimensional disordered systems studied till now these two 
quantities always coincide). 

Both M- and H-stacks demonstrate a rather narrow 
crossover from the localized to the ballistic regime. The 
H-stack in the near ballistic region, and the M-stack 
in the ballistic region are weakly scattering disordered 
stacks, while in the far ballistic region, the H-stack trans- 
mits radiation as an effectively uniform medium. 

We also consider the effects of loss and show that ab- 
sorption dominates the effects of disorder at very short 
and very long wavelengths. The crossover region is par- 
ticularly sensitive to losses, so that even small absorp- 
tion suppresses oscillations in the transmission length as 
a function of frequency. 

All of the theoretical results mentioned above are con- 
firmed by, and are shown to be in excellent agreement 
with, the results of extensive numerical simulation. 

In M-stacks with only refractive index disorder, An- 
derson localization and transmission resonances are ef- 
fectively suppressed and the crossover region between 
the localized and ballistic regimes is orders of magnitude 
greater. A more detailed study of the A^-anomaly shows 
that the genuine wavelength dependence of the transmis- 
sion length is not described by any power law and rather 
is non-analytic in nature. 

In what follows. Sec. HIl presents a detailed description 
of our model. Section IIIII is devoted to the analytical 
studies of the problem, while the results of numerical 
simulations and a discussion of these are presented in 
Sec.IVl 



II. MODEL 

A. Mixed and homogeneous stacks 

We consider a one-dimensional alternating M-stack, as 
shown in Fig. [TJ It comprises disordered mixed left- (L) 
and right- (R) handed layers, which alternate over its 
length of N layers, where N is an even number. The 
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thicknesses of each layer are independent random vahies 
with the same mean value d. In what follows, all quanti- 
ties with the dimension of length are measured in units 
of d. In these units, for the thicknesses of a layer we can 
write 
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Although localization in disordered H-stacks with 
right-handed layers has been studied by many 
authoraiSii^i^^""— , here we consider this problem in 
its most general form and show that the transmission 
properties of disordered H-stacks are qualitatively 
the same for stacks comprised of either solely left- or 
right-handed layers. 



where the fluctuations of the thickness, S'^'^\ j 
are zero-mean independent random numbers. 
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FIG. 1: (Color online) Structure geometry. 

the magnetic permeability for right-handed media to be 
fij — 1 and for metamaterials to be = — 1, while the 
dielectric permittivity is 



e,=±(l + 4'')±^a,)^ 

where the upper and lower signs respectively correspond 
to normal (right-handed) and metamaterial (left-handed) 
layers. The refractive index of each layer is then 



B. Transmission: localized and ballistic regimes 

We introduce the transmission length In of a finite 
random configuration as 



1 

In 



In IT, 



N\ 



N 



(1) 



where T/v is the random transmission coefficient of a sam- 
ple of the length N. As a consequence of the self-averaging 
of ln|T^|/7V, 



lim 



InlT^I ,. 1 1 
— lim — = -. 

N N^oo In I 



(2) 



This means that for a sufficiently long stack (in the local- 
ized regime) the transmission coefficient is exponentially 
small iTjvl ~ exp{-N/l). 

In what follows, we consider stacks composed of lay- 
ers with low dielectric contrasts, i.e, IS"^'"^] <C 1, so that 
the Fresnel reflection coefficients of each interface, and of 
each layer, are much smaller than 1. Here, a thin stack 
comprising a small number of layers is almost transpar- 
ent. In this, the ballistic regime, the transmission length 
takes the form 



where all jj^'* and absorption coefficients of the slabs, 
(Tj > 0, are independent random variables. With this, 
the impedance of each layer relative to the background 
(free space) is 



with the same choice of the sign. 

We begin with the general case when both types of dis- 
order (in refractive index and in thickness) are present. 
Two particular cases, each with only one type of disor- 
der, are rather different. In the absence of absorption, 
the M-stack with only thickness disorder is completely 
transparent, a consequence of Zj = 1. However, the case 
of only refractive-index disorder is intriguing because, as 
is shown below, such mixed stacks manifest a dramatic 
suppression of Anderson localization in the long wave 

2? 

regiortsi. 



1_ ^ {\Rn\') 
In^ 2N ' 

involving the average reflection coefficient^, which is 
valid in the case of lossless structures. This follows di- 
rectly from Eq. ([1} by virtue of the conservation relation- 
ship, |i?Arp -I- |T/vp = 1. Thus, in the ballistic regime, 

9N 

, iV«6. (4) 

where the length b in this equation is termed the ballistic 
length. 

Accordingly, in studies of the transport of the classical 
waves in one-dimensional random systems, the following 
spatial scales arise in a natural way: 

• In — the transmission length of a finite sample ((ij , 

• I — the localization length and 

• b — baUistic length (g]). 
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Note that in the case of absorbing stacks {aj = ct > 0), 
the right hand side of Eq. ([T]) defines the attenuation 
length, latt , which incorporates the effects of both disor- 
der and absorption. 

In what follows, we show that contrary to commonly 
accepted belief, the quantities 7~^, /, and b are not nec- 
essarily equal, and, under certain situations, can differ 
noticeably from each other. 

In this paper, we study mainly the transmission length 
defined above by Eq. [T] This quantity is very sensi- 
tive to the size of the system and therefore is best suited 
to the description of the transmission properties in both 
the localized and ballistic regimes. More precisely, the 
transmission length coincides either with the localization 
length or with the ballistic length, respectively in the 
cases of comparatively thick (localized regime) or com- 
paratively thin (ballistic regime) stacks. That is. 



Here, /J, 



k = 2tt/X, and A denotes the dimen- 



b N <$:b. 



Another argument supporting our choice of the transmis- 
sion length as the subject of investigation is that it can 
be found directly by standard transmission experiments, 
while measurements of the Lyapunov exponent call for a 
much more sophisticated arrangement. 



sionless free space wavelength. While the sign of the 
phase shift across each slab Re {f3j ) varies according to 
the handedness of the material, the Fresnel interface co- 
efficient pj given by 



Pj = 



Z, + 1 ' 



(9) 



depends only on the relative impedance of the layer Zj , 
a quantity whose real part is positive, irrespective of the 
handedness of the material. 

Equations ©-([HI) are general and provide an exact de- 
scription of the system and will be used later for direct 
numerical simulations of its transmission properties. 

It was mentioned previously that we consider the spe- 
cial case of weak scattering for which the reflection from a 
single layer is small, i.e., \rj \ ^ 1. This occurs either for 
weak disorder, or for strong disorder provided that the 
wavelength is sufficiently long. The transmission length 
then follows from 



IuITatP = 2RelnrAr, (10) 

and requires the following first order approximations de- 
rived from Eqs. ([5]) and ([6]): 



III. ANALYTICAL STUDIES 



A. Weak scattering approximation 



lnT„ = lnri,„_i + lnt„ + i?„_ir„, (11) 
Rn = r„ + Rn^itl. (12) 



The theoretical analysis involves the calculation of 
the transmission coefficient using a recursive procedure. 
Consider a stack which is a sequence of N layers enu- 
merated by index n from n = 1 at the rear of the stack 
through to rt = at the front. The total transmis- 
sion (r„) and refiection (i?„) amplitude coefficients of 
the stack satisfy the recurrence relations 



Tn = 
Rn = 



Tfi — ltn 

1 - Rn-lTn ' 

n 4-2 



1 - Rn-lTn 



(5) 
(6) 



for n — 2, . . . ,N, in which both the input and output 
media are free space. In Eqs ^ and ©, the amplitude 
transmission (tj) and the refiection (rj) coefficients of a 
single layer are given by 



Pj-(l-e''ft) 



t, = 



l-p2g2z,3, ' 

(1 - P,^)e'^^ 
1 - p2e2^ft ■ 



(7) 

(8) 



In deriving Eq. ()12p . we omit the first-order term 
Rn-i^n^n since it contributes only to the second order 
of In Tn already after the first iteration. Then, by sum- 
ming up logarithmic terms (jlip . we obtain 



TV 



N N 



InT, 



TV 



n 



tl (13) 



m— 2 j—m 



p—j—m+2 



This equation enables us to derive a general expres- 
sion for the transmission length It{N) which is valid in 
all regimes (see the next Section). However, the ballistic 
length b, according to Eq. ^ (and the average reflec- 
tion coefficient as well), is determined only by the total 
reflection amplitude Rn in the case of lossless strutures. 
In the ballistic regime, this amplitude coefficient, to the 
necessary accuracy, is given by 



R 



N 



N 

E^ 



(14) 
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B. Mixed stack 

1. General Approach 

Prom here on, we assume that the random variables 
5'^^'' of left-handed or right-handed layers, and ctj 
are identically distributed according to the correspond- 
ing probability density functions. This enables us to ex- 
press all of the required quantities via the transmission 
and reflection amplitudes of a single right-handed or left- 
handed layer, tr,u ^r,i, and also to calculate easily all of 
the necessary ensemble averages. 

The average of the first term in Eq. can be written 

as 



Ini,) = -{lntr) + -{\nti). 
Next, we split the second term of Eq. into two parts 

N N j-1 

51 '^j-'^+i'^j n tl = Nni + Nn2, 

m—2j—7n p—j—m+2 

where 



In the lossless case (cr = 0), the first term on the right 
hand side of Eq. corresponds to the so-called single- 
scattering approximation, which implies that multi-pass 
reflections are neglected so that the total transmission 
coefficient is approximated by the product of the single 
layer transmission coefficients, i.e., 

N 

In the case of very long stacks (i.e., as the length 
N — > 00), we can replace the arithmetic mean, 

Y^^=i In \tj\, by its ensemble average (In |t|). On the 
other hand, in this limit the reciprocal of the transmission 
length coincides with the localization length. Using the 
energy conservation law, |rjp -I- jt^ p = 1, which applies 
in the absence of absorption, the inverse single-scattering 
localization length may be written as 

and is proportional to the mean reflection coefficient of a 
single random laye r-'^'^'^" . The corresponding modification 
of Eq. HH) then reads 



7^l 



7^2 = 



^ N/2 N 



ti 



i-1 

m— 1 j—2m p—j — 2ni-\-2 

^ N/2 N j-1 

m— Ij— 2?n — 1 p— j — 2m+l 



comprising contributions to the depletion of the transmit- 
ted field due to two pass reflections respectively between 
slabs of different materials (i.e., of opposite handedness), 
and between slabs of the same material (i.e., of like hand- 
edness). Averaging these expressions, we obtain 



where 



= ArAi 



1 



{B^ -1)(1 + S2) 



52 



2N 



iV(l-B2)2 

N 2{B^ - 1 



1-52 ' (1-52)2 



, (15) 
(16) 
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In 



{\ri\ 



Re((7^l) + (7^2)). 



(20) 



Here, the first term corresponds to the single-scattering 
approximation, while the next two terms take into ac- 
count the interference of multiply scattered waves as well 
as the dependence of the transmission length on the stack 
size. Note that Eq. ([20)1 is appropriate only for lossless 
structures. In the presence of absorption, Eq. ([T8| should 
be used instead. 



2. Transmission length 

From this point on, we assume that the statistical prop- 
erties of the right-handed and left-handed layers are iden- 
tical. As a consequence of this symmetry, the following 
relations hold for any real- valued function g in cither the 
lossless or absorbing cases: 



Ar^{rr), Br^{tl), T^l,r, B^^BiB,.. (17) 

The resulting transmission length is determined by the 
equation 

-i-= "°l'-l':"°l"l' H-Re((K.)H-(7^.». (18) 
In ^ 



{giU)) - {g{U)r, {gir,.)) = {gin))*. (21) 

Therefore, Ar and Br arc the complex conjugates of Ai 
and Bi, and B^ is real quantity, as are both averages 
(7^l) and (7^2). 

Accordingly, as a consequence of the left-right symme- 
try (pij) . the transmission length of a M-stack depends 
only on the properties of a single right-handed layer and 
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may be expressed in terms of three averaged characteris- 
tics: (r), (hi and (t^) (in which we omit the subscript 

r)- ^ 

With these observations, the transmission length of a 
finite length M-stack may be cast in the form: 



1 

In 



fiNj). 



(22) 



where 



1 . ... \{r}\^+Re{{r)^{tr) 
7 = -aH^I> Y^J^, : 



(23) 



and 



2/1 



l-exp(-2/0 
'\{r)\^+Re{{rf{f 



\{r)\' 



(24) 



are, as we will see below, the inverse localization and 
inverse ballistic lengths. The function f{N, I) is defined 
as 



f{Nj) 



I 

N 



exp 



N 

T 



(25) 



and introduces a new characteristic length termed the 
crossover length 



(26) 



In KOI' 



which arises in the calculations in a natural way and, 
as will be demonstrated below, plays an important role 
in the theory of the transport and localization in one- 
dimensional random systems. Equations (|22l) - (p4|) com- 
pletely describe the transmission length of a mixed stack 
in the weak scattering approximation. 

Obviously, the characteristic lengths Z(A), b{X), and 
1{X) appearing in Eq. ([2^ are functions of wavelength. 
Using straightforward calculations, it may be shown that 
the first two always satisfy the inequality 1{X) > b{X), 
while, as we will see, the crossover length is the shortest 
of the three, i.e., b{\) > /(A) in the long wavelength 
region. 

In the case of a fixed wavelength A and a stack so short 
that N ^ '(A), the expansion of the exponent in Eq. 
p5|) yields / — ^ 1, in which case the transmission length 
approaches b{X). Correspondingly, for a sufficiently long 
stack, N ^ '(A), / and the transmission length 
assumes the value of 1{X). 



In summary. 



InW 



1{X), A^>Z"(A), 
6(A), 7V<r(A), 



with the transition between the two ranges of iV being 
determined by the crossover length. 

While in the lossless case, the ballistic regime occurs 
when the stack is much shorter than the crossover length 
(iV <C /(A)), it is important to note that, in the local- 
ization regime, the opposite inequality is not sufficient 
and the necessary condition for localization is N ^ Z(A). 
In what follows, we consider samples of an intermediate 
length, i.e., 1{X) < iV < Z(A). 

For a M-stack of fixed size N, the parameter governing 
the transmission is the wavelength, and the conditions for 
the localized and ballistic regimes should be formulated 
in the wavelength domain. To do this, we introduce two 
characteristic wavelengths, Ai(A^) and X2{N), defined by 
the relations 



N = liX,{N)), N^l{X2{N)). 



(27) 



It can be shown that the long wavelength region, A ^ 
Ai(A^), corresponds to localization where the transmis- 
sion length coincides with the localization length, while 
in the extremely long wavelength region, A 3> A2(A'^), 
the propagation is ballistic, with the transmission length 
given by the ballistic length b. That is. 



In{X) 



1{X), A«Ai(iV), 
6(A), A»A2(iV). 



(28) 



When Ai (N) < X2 (N), there exists an intermediate range 
of wavelengths, Ai(A^) < A < X2{N), which will be dis- 
cussed below. 

To better understand the physical meaning of the ex- 
pressions and for the localization and ballis- 
tic lengths, we will consider ensembles of random con- 



figurations in which the fiuctuations S^'^'^ and S'^"'' are 
distributed uniformly over the intervals [—Q^,Qu] and 
[—Qd, Qd] respectively, with aj = 0. The average quanti- 
ties that arise in Eqs. ((23|) and (l24l) are presented in the 
Appendix. These formulae allow for calculations with an 
accuracy of order 0{Q^) for arbitrary Q^. For the sake 
of simplicity, we assume that the fluctuations of the re- 
fractive index and thickness are of the same order, i.e., 
Qv ^ Qd so that the dimensionless parameter 
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Qi 



is of order of unity. We also neglect the contribution of 
terms of order higher than Q?. 



7 



The short wavelength asymptotic of the locaUzation 
length is then 



12 



(29) 



In the long wavelength limit, we obtain the following 
asymptotic forms for the corresponding single layer aver- 
ages: 



(ln|t|) 



6 2 



6 ' 



1 + 2ifc 



ikQl 



-2k^ - 



(31) 



(32) 



Substitution of these expansions into Eq. (|24p yields the 
long wavelength asymptotic of the ballistic length 



6(A) 



3A2 



(33) 



This asymptotic can be calculated directly from Eq. 
dH) with the average reflection coefficient, determined 
from Eq. p^ . being 



\RNf)=N[{\rf)-\{r)f)+N 



2/Rer)^ 



(34) 



The same substitutions into this equation (|34|) give the 
average total reflection coefficient 



R 



N\ 



Nk^Ql 



(35) 



In the ballistic regime, the final term is negligibly small. 
This, together with Eq. again results in the value of 
the baUistic length given in Eq. 

Substituting the long wavelength expansions (l30l) - (|32)) 
into Eqs. and (pS)) . we derive the following asymp- 
totic forms for the localization length 



3A2 3 + C 



(36) 



and the crossover length 



3A^ 



27r2g2 4(3 + 0' 



(37) 



In seeking to compare the result (136)) for the localiza- 
tion length with the corresponding long wave asymptotic 
of the reciprocal Lyapunov exponent, it is important to 
note that these two quantities are defined in different 



ways. The localization length is defined by Eq. ([2]) via 
the transmittivity on a realization, while the Lyapunov 
exponent describes an exponential growth of the enve- 
lope of a currentless solution far from a given point in 
which the solution has a given valuei^. We have calcu- 
lated the long wave asymptotic of the M-stack Lyapunov 
exponent using the well known transfer matrix approach, 
assuming the same statistical properties of the modulus 
of dielectric constant and the thickness of both left- and 
right-handed layers. 



■K^d? e2 - £2 
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2A2 e ' 



(38) 



In the case of rectangular distributions of the fiuctuations 
of the dielectric constants and thicknesses, this result re- 
duces to 
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27r2Q2 
3A2 



1 3 + C 1 
/ 1 + C ^ ^ 



(39) 



(we neglected the small corrections proportional to <C 
1). Thus, the disordered M-stack in the long wavelength 
region presents a unique example of a one-dimensional 
disordered system in which the localization length differs 
from the reciprocal of the Lyapunov exponent. 

The long wavelength asymptotic of the Lyapunov ex- 
ponent 7 (X A^2 aiifj the asymptotic of the localization 
length, Z (X A2 , Eq. ([55]) have rather clear physical mean- 
ing. Indeed, in the limit A — >■ cx), the propagating wave 
is insensitive to disorder since the disorder is effectively 
averaged over distances of the order of the wavelength. 
This means that Anderson localization is absent and the 
Lyapunov exponent 7 cx vanishes. For large but fi- 
nite wavelengths (i.e., small wavenumbers k = 27r/A), 
the Lyapunov exponent is small and, assuming that its 
dependence on the wavenumber is analytic, we can ex- 
pand it in powers of k. This expansion commences with a 
term of order k"^ since the Lyapunov exponent is real and 
the wavenumber enters the field equations in the form 
(ifc). Accordingly, in the long wavelength limit, 7 cx fc2 
and so / cx A2. 

The behavior of the Lyapunov exponent given by (j38p 
does not depend on the handness of the layers; compare 
this equation with the corresponding long wave asymp- 
totic obtained in Refs.^^'25 foj- homogeneous stacks com- 
posed of only positive or only negative layers. However, 
it crucially depends on the propagating character of the 
field within a given wavelength region. When this re- 
gion is within the gap of the propagation spectrum of 
the effective ordered medium, the frequency dependence 
of the Lyapunov exponent differs from ^ A~2. It takes 
place, for example, for a stack composed of single nega- 
tive layers, where only one of two characteristics (e, /i) is 
negative25.. 

The main contribution to the ballistic length is 
due to the final term in the right hand side of Eq. 
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which corresponds to the single-scattering approximation 
discussed at the end of the Sec. IIIIB II While the bal- 
listic length follows from the single scattering approxi- 
mation, the calculation of the localization length (1551) is 
more complex and requires that interference due to the 
multiple scattering of waves must be taken into account. 

In the case under consideration (Q^, ^ Qd and C ^ 1), 
the two characteristic wavelengths Ai(A^) and A2(A^) ((27)) 
for an M-stack of a fixed length size N take the form 



Ai(iV) = ttQ,. 



/2A^14-C 
~3 + C' 



X2iN) = 27tQ^. 2N 1 



3/' 



(40) 
(41) 



and are of the same order of magnitude. This means that 
the transmission length In coincides with the localization 
length (|36|) for N ^ '(A)), and with the ballistic length 
([55| in the case when N < [(A). 

The crossover between the localized and ballistic 
regimes occurs for Ai(A^) ;$ A < A2(A^), with the two 
bounds being proportional to Q^\/N and thus growing 
as ^/N. For long stacks, the transmission length in the 
crossover region can be described with high accuracy 
by the general equations (|22|) and (|25|) . where the bal- 
listic length 6(A), the localization length /(A), and the 
crossover length ^(A) are replaced by their asymptotic 
forms (I55)) - (|T7|) . In the case of solely refractive index 
disorder (Qd 0) when all thicknesses are set to unity, 
the ballistic length and the short wavelength asymptotic 
of the localization length coincide with the limiting values 
of their counterparts for an M-stack with both refractive 
index and thickness disorder ([33]) and (j29p . However, the 
corresponding limiting values of the long wavelength lo- 
calization length and crossover lengths have nothing to 
do with the genuine behaviour of the transmission length 
(see Sec. IIVB II) . This means that the weak scatter- 
ing approximation fails to describe the long wavelength 
asymptotics of the transmission length in both the local- 
ization and crossover regions. This is discussed in greater 
detail below. 

Eqs. (dHj-dlHl) (or Eqs. (HHj-dH])) in the symmetric 
case) completely determine the behaviour of the trans- 
mission length for a mixed stack composed of weakly 
scattering layers. Although I has been introduced as a 
crossover length, the entire region N does not neces- 
sarily support localization. Correspondingly, the ballistic 
regime may exist outside the region N <^ I. 

All of the results obtained under the symmetry as- 
sumption (j2ip are qualitatively valid in the general case. 
Indeed, the existence of the crossover (P5| is related to 
the exponential dependence in Eq. (j25|) with a real and 
positive crossover length I. When the assumption of sym- 
metry no longer holds, this length takes a complex value. 
However, the quantity B in Eq. (IT51) satisfies (by its 
definition in Eq. (|17l) ) the evident inequality \B\ < 1. 



Therefore, the corresponding analogue of the function / 
(PSj) preserves all necessary limiting properties (see the 
next Sec. lIHBSI) . 



3. Homogeneous stacks 

In this Section, we consider a H-stack composed en- 
tirely of normal material layers noting that the behaviour 
of a H-stack of metamaterial (left-handed) layers alone is 
exactly the same, a result which may be obtained directly 
from Eq. (jl8l) by replacing each I by r, after which any 
reference to the index r may be omitted. The transmis- 
sion length of an H-stack is then 



1 

In 



1 

N 



Re 



(r) 



2\N 



1 - it') 
(l-(t2))2 



where the inverse localization length I is 



j = -(ln|i|: 



Re 



I -it')' 



(42) 



(43) 



Now we consider a H-stack composed of weakly scatter- 
ing layers. To simplify the discussion, we consider only 
refractive index disorder (i.e., Qd — 0). In this case, 
the asymptotic behavior of the localization length in the 
short and long wavelength limits are 



^2 A 



Z(A) 



3A2 



(44) 



A 



The main contribution to the localization length is re- 
lated to the first term in Eq. (|43| . Thus, the localization 
length of the H-stack in the long wavelength region is 
described completely by the single scattering approxima- 
tion and coincides with the ballistic length p3p of the 
M-stack. 

Using the transfer matrix approach, we can also cal- 
culate the long wavelength asymptotic of the Lyapunov 
exponent for a H-stack . It is described by the same equa- 
tion ([55]) as the asymptotic for the M-stack, thus coincid- 
ing with the asymptotic of the reciprocal Lyapunov expo- 
nent. In recent work22,, this coincidence was established 
analytically in a wider spectral region. However, the nu- 
merical calculations intended to confirm this result are 
rather unconvincing. Indeed, the numerically obtained 
plots demonstrate strong fluctuations (of the same order 
as the mean value) of the calculated quantity, while the 
genuine Lyapunov exponent is non-random and should 
be smooth without any additional ensemble averaging 
mentioned by the authors. 

If we cast (t^)^ in the second term of Eq. (^^ in the 
form exp (A^ln(i^)) we see that the crossover length of 
the H-stack is 
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1= |ln(t2)| \ 

with its long wavelength asymptotic, according to Eq. 
(IMl), being 



[(A) 



X_ 

An' 



(45) 



Here, the crossover length is proportional to the wave- 
length, in stark contradistinction to the situation for M- 
stacks, in which the crossover length is proportional to 

To consider this further, we define the characteristic 
wavelengths Ai(A^) and X2{N) by the expressions (l?7l) . 
For a H-stack, these lengths are 



\2{N) = AnN. 



(46) 



Evidently, the second characteristic wavelength is always 
much larger than the first X2{N) ^ Ai(A^). As a con- 
sequence, the long wavelength region, where the ballis- 
tic regime is realized, can be divided into two subre- 
gions. The near subregion (moderately long wavelengths) 
is bounded by the characteristic wavelengths 



Ai(iV) < A < A2(iV). 

The main contribution to the ballistic length in the near 
ballistic region, 6„, is due to the first term in Eq. (j43| . 
Thus the ballistic length &„(A) has the same wavelength 
dependence as the localization length 1{X) (|44l) 



2n4 



bfW 



3A2 



18A2 



(47) 



In the case of a relatively short H-stack NQf,/ 12 ^ 1, 
the contribution of the first term in the right hand side of 
this equation dominates, and hence the transition from 
the near to the far subregions is not accompanied by any 
change in the analytical dependence on the wavelength. 
The ballistic length is thus described by the same wave- 
length dependence over the entire ballistic region 



&(A) = 



3A2 



27r2g2 



Ai(A^) < A. 



(48) 



For sufficiently long H-stacks NQ^/ 12 ^ 1, in the far 
long wavelength region, the second term is dominant and 
so 



18A2 



(49) 



Thus, the wavelength dependence of the ballistic length 
of a sufficiently long H-stack is 



6(A) 



3A2 



18A2 



Ai(iV)<A< A2(iV), 



>^2{N)<X. 



(50) 



The same result for the ballistic length also follows 
from Eq. ^ with Eq. ([HI) , in the case of a homogeneous 
stack, yielding 



bniX) 



3A2 



and is well described by the single scattering approxima- 
tion ([T^ . For a H-stack , the transition from the local- 
ized to the ballistic regime at A ~ Ai(A^) is not accom- 
panied by any change of the wavelength dependence of 
the transmission length. This change can occur at much 
longer wavelengths A '-^ A2 (N) in the far long wavelength 
subregion. 

To derive the ballistic length 6/ in the far long wave- 
length region, we may proceed in a similar manner to 
that outlined in the case of a M-stack and expand the 
exponent (t^)^ = exp (A^ln(t^)) in Eq. gH). However, 
because (i^) in this expression is complex, the situation 
is more complicated than was the case for the M-stack. 
In particular, the first two terms of the expansion do 
not contribute to the ballistic length. Taking account of 
the second order leads to the following expression for the 
ballistic length, 6/ (A), in the far long wavelength region 



[\RN\') = N[(\r\')-\{r)\')+N'\{r)\\ (51) 
In the long wave limit this leads to 



\Rn\' 



3 36 ^ 



Ai(7V)«A, (52) 



The final term in Eq. (1511) differs from the final term 
in Eq. (|34l) and, therefore, in contrast to the M-stack, its 
contribution to the total reflection coefficient can be of 
the order of, or larger than, that of the flrst term. This 
together with Eq. (jlj is equivalent to the result in Eqs. 
(|48l) and ([50l) . and is applicable to short and long stacks 
respectively. 

The far long wavelength ballistic asymptotic (gS]) has a 
simple physical interpretation. Indeed, in this subregion, 
the wavelength essentially exceeds the stack size and so 
we may consider the stack as a single weakly scattering 
uniform layer with an effective dielectric permittivity ScB- 
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In this case, the bahistic length of the stack according to 
Eq.dll) is 



bf 



2N 



toff 



(53) 



where 



-Rcff 



ikN 



(Eoff - 1) 



(54) 



is the long wavelength form of the reflection amplitude 
for a uniform right-handed stack of the length N and 
constant dielectric permittivity ScS (with fj, — 1). In 
the case of a uniform left-handed stack, the value of the 
reflection coefficient should be 



RcS 



ikN 



(ecff -f- 1) : 



(55) 



Eq. ([58]), together with Eqs ([53]), ([54]), leads imme- 
diately to the far long wavelength ballistic length (|49p . 
We emphasize that because of the effective uniformity 
of the H-stack in the far ballistic region, the transmis- 
sion length on a single realization is a less fluctuating 
quantity. In contrast, transmission length for H-stacks 
fluctuates strongly in the near ballistic region, as indeed 
it does over the entire ballistic region for M-stacks. 

To characterize the entire crossover region between the 
two ballistic regimes ([5(7)) in greater detail, we return to 

the general formula ([1^ . in which we represent (t^)^ as 
exp(iV In (t^)). Then, in accordance with Eq. ([5^ . we 
can write 

(t^)^ « exp{N{2ik - k^Ql) « (1 - Nk^Ql) exp{2ikN). 

(61) 

As a consequence, instead of the result in Eq. ([47|. we 
obtain 



in which the value of £off is now negative (with fi = — 1)- 
To calculate the effective parameters, we use Eq.([T3]). 
Neglecting small fluctuations of the layer thickness, the 
single layer reflection amplitude in the long wavelength 
limit reads 



rj = y {ej - 1) ■ 



(56) 



The total reflection amplitude for a stack of length N is 
ik , , ikN I 1 1 

^^ = yEfe-i) = ^hvE^^-i ' (57) 



corresponding to the effective dielectric permittivity de- 
termined by the expression 



1 ^ 



(58) 



i=i 



For sufficiently long stacks, the right hand side of this 
equation can be replaced by the ensemble average {ej) 
and hence 



Ecff — (£7 



(59) 



This expression is similar to that of the two-dimensional 
case for disordered photonic crystals reported in Ref.^^. 
The same considerations as above, but for a homoge- 
neous stack composed entirely with disordered metama- 
terial slabs, also lead to an effective value of the permit- 
tivity. 



Eeff 



(60) 



&n(A) 



3A2 



727V 



27rN 
A ' 



(62) 



We note that, strictly speaking, the expansion ([6T]) is 
vahd inside the interval V6Xi{N) < A < A2(iV). 

The second term in Eq. ([5^ represents standard os- 
cillations of the reflection coefficient of a uniform slab of 
finite size. In the far long wavelength limit A ^ A2, this 
equation coincides with (|47p . Thus, Eq. ([5^ describes 
the ballistic length over practically the whole ballistic re- 
gion A ^ Ai(iV). Moreover, taking into account that the 
long wavelength asymptotic of the localization length co- 
incides with that of the near long wavelength ballistic 
length, we see that the right hand side of Eq. ([52]) serves 
as an excellent interpolation formula for the reciprocal 
of the transmission length l/Zjv(A) of a sufficiently long 
stack {NQf, /12 ^ 1) over the entire long wavelength re- 
gion A > 1. 



C. Comparison of the transmission length 
behaviour in M- and H- stacks 

Away from the transition regions, the transmission 
length can exhibit three types of long wavelength asymp- 
totics described by the right hand sides of Eqs. ([55]) , ([5S]) , 
and (|49p . The flrst ([55]) corresponds to the single scatter- 
ing approximation where the inverse transmission length 
is proportional to the average reflection coefficient of a 
single random layer. In the absence of absorption, this 
characterizes both the localization and ballistic lengths. 
The second asymptotic form (|36p takes into account the 
interference of multiply scattered waves and describes 
the localization length. The third asymptotic (H^ cor- 
responds to transmission through a uniform slab with 
an effective dielectric constant given by Eq. ([59| and is 
relevant only in the ballistic regime. 

In the case of a M-stack, the first two expressions ([55]) 
and (1551) characterize the ballistic and localized regimes 
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respectively, while the third asymptotic is never realized 
in M-stacks. In relatively short H-stacks the long 

wavelength behaviour of the transmission length is de- 
scribed by the same dependence (1551) in both the local- 
ized and ballistic regimes. Finally, in the case of long 
H-stacks (j49|). the transmission length follows from Eq. 
P3p in the localized and near ballistic regimes, while in 
the far ballistic region it is described by the right hand 
side of Eq. p9|) . 

These results predict the existence of a different wave- 
length dependence of the transmission length in differ- 
ent wavelength ranges. For M-stacks, the crossover be- 
tween them occurs at the wavelength Ai(A^) (l40t where 
the size of the stack is comparable with the localization 
length l{Xi{N)) « N. For long H-stacks, the crossover 
occurs when the wavelength becomes comparable to the 
stack size X ^ N. Short H-stacks exhibit the same wave- 
length dependence of the transmission length in all long 
wavelength regions, i.e., in both localized and ballistic 
regimes. 
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FIG. 2: (Color online) Transmission length Zat vs A for 
M-stack (thick solid line) and H-stack (thick dashed line). 
Asymptotics of the localization length I (thin straight lines), 
the short wavelength asymptotic (thin dotted line, ((29]), and 
the long wavelength asymptotics — thin sohd hne for the M-- 
stack (|36p and a thin dashed line for the H-stack (|44p . 



IV. NUMERICAL RESULTS 

The results of the numerical simulations presented be- 
low correspond to uniform distributions of the fluctua- 
tions J^'^^ 5^'"^ with widths of Q^, and Qd, respectively, 
and with a constant value of the absorption a in each 
layer. 

Results are presented for (a) direct simulations based 
on the exact recurrence relations ([5]) and ([6]); (b) the weak 
scattering analysis for the transmission length based on 
Eqs. dm and (jtH); and (c) asymptotics for short and 
long wavelengths [ (P5)) . (I55|) . and (|^ ] . For the mean 
reflectivity, we used the asymptotic forms (|35|) and (|52]) . 

In all cases, unless otherwise is mentioned, the ensem- 
ble averaging is taken over Nj. = realizations. The 
results up to Sec. IIVCI are for lossless stacks (a = 0) 
only. 



A. Refractive-index and thickness disorder 

We first consider stacks having refractive index and 
thickness disorder, with = 0.25 and Qd = 0.2. Shown 
in Fig. [5] are transmission spectra for a M-stack of = 
10^ layers and a H-stack of length N = 10'^. There are 
two major differences between the results for these two 
types of samples: first, in the localized regime {N ^ 
In), the transmission length of the M-stack exceeds or 
coincides with that of the H-stack; second, in the long 
wavelength region, the plot of the transmission length of 
the M-stack exhibits a pronounced bend, or kink, in the 
interval A G [10^, 10'^], while there is no such feature in 
the H-stack results. These two types of behaviour are 
discussed in more detail below. 



1. Mixed stacks 

The weak scattering approximation (WSA) of Eq. (fT5)) 
is an excellent method by which to calculate the trans- 
mission length for M-stacks. This is seen in Fig. [2] where 
the curves obtained by numerical simulations and by the 
WSA are indistinguishable (solid line). The character- 
istic wavelengths (|40)) of this mixed stack are Ai ~ 148 
and A2 ~ 839. Therefore, the transmission length de- 
scribes the localization properties of a random sample 
in the region A < 148, whereas longer wavelengths, A > 
839, correspond to the ballistic regime. The crossover 
from the localized to the ballistic regime demonstrates 
the kink-type behaviour that occurs within the region 
•^1 ^ ^ ^2- The short and long wavelength behaviour 
of the transmission length is also in excellent agreement 
with the calculated asymptotics in both regimes. 

To analyze the long wavelength region (A > 10) more 
carefully, we plot in Fig.[3]the transmission lengths of M- 
stacks of three different sizes, N — 10"^, 10^, and 10^. In 
all cases, there is excellent agreement between the sim- 
ulations and the WSA predictions. The characteristic 
wavelengths for TV = 10^ are Ai = 14.8 and A2 = 83.9, 
while for N = lO'' they are Ai = 1480 and A2 = 8390, and 
we see that the observed crossover regions are bounded 
exactly by these characteristic wavelengths in all three 
cases. 

To confirm the ballistic nature of the transmission in 
the region A > A2, we plot in Fig. HJa) the logarithm 
of the mean value of the reflectance for the same three 
stack sizes as a function of the logarithm of the wave- 
length. In all cases, the plots exhibit a linear dependence 
ln(|i?Arp) = const -|- 21nA in the ballistic regime, which 
is bounded from below by the crossover wavelength A2. 
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FIG. 3: (Color online) Transmission length In vs \ for a M- 
stack of A'^ = 10^ (thick solid line), 10^ (thick dashed line) and 
lO'^ (thick dotted line) layers showing both numerical simula- 
tions and the WSA theory. The long wave asymptotics for the 
localization length p6p and the ballistic length (|33p are shown 
respectively in the thin solid and dashed lines respectively. 



The straight lines are calculated from Eq. ([35|) and con- 
firm that the refiection coefficient in the mixed stack is 
proportional to the stack length. Within the localized 
region A < Ai, the reflection coefficient is close to unity 
in all three cases. 

The behaviour of the transmission length is illustrated 
by the phase diagram in the (A, N) plane shown in 
Fig. Iljb). The two slanted fines N ^ 1{X) and N = 1{X) 
separate the plane into three parts corresponding to the 
localization (I) , the crossover region (II) and the ballistic 
region (III) . The intersections of these lines with the hor- 
izontal lines N ^10^, N = 10^, and N = 10^ define the 
characteristic wavelengths Ai and A2 for the three stack 
sizes considered here. It is easy to see that these wave- 
lengths, determined with the aid of the phase diagram, 
perfectly bound the crossover regions in both Fig. |3]and 
Fig. Ha). 

Until now, we have dealt only with the transmission 
length /at (A), which was defined through an average 
value. However, more detailed information can be ob- 
tained from the transmission length ?jv(A) for a single 
realization, defined by the equation 
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FIG. 4: (Color online) (a) Average reflectance for M-stacks 
of length TV = 10^ (solid line), lO'^ (dashed line) and 10^ (dot- 
ted line) layers (numerical simulation and WSA). Long wave 
asymptotic for the average reflectance for the same stacks 
(thin solid lines). 

(b) Phase diagram of M-stacks. The thick solid line corre- 
sponds to a stack size equal to the localization length. The 
dashed line corresponds to a stack size equal to the crossover 
length. The localized, crossover and ballistic regimes occur in 
regions I, II, and III respectively. 



1 ^ In iTjvl 



(63) 



In the localized regime, for a sufficiently long (i.e., 
N ^ I) M-stack, the transmission length for a single re- 
alization In{X) is practically non-random and coincides 
with In{X), while in the ballistic region it fiuctuates. The 
data displayed in Fig. [S]) enables us to estimate the dif- 
ference between the transmission length In{X) (solid line) 
and the transmission length InW for a single randomly 
chosen realization (dashed line), and the scale of the cor- 
responding fluctuations. Both curves are smooth, coin- 
cide in the localized region, and differ noticeably in the 
ballistic regime. The separate discrete points in Fig. [S] 
present the values of the transmission length Zjv(A) cal- 
culated for different randomly chosen realizations. It is 
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FIG. 5: (Color online) Transmission lengths In (solid black 
line) and the transmission length for a single realization Im 
(dashed blue line) vs A for a M-stack with Qi, — 0.25, Qd ~ 0.2 
and = 10^ layers. Each separate point corresponds to a 
particular wavelength with its own realization of a random 
stack. 



10" 



evident that fluctuations in the ballistic region become 
more pronounced with increasing wavelength. 



2. Homogeneous stacks 

The absence of any kink in the H-stack transmission 
length spectrum in Fig. [2] follows from Sec. IIII B 31 in 
which it was shown that the crossover to the far bal- 
listic regime occurs at the wavelength X2{N) pS)) . For 
N = lO'^, this is of the order of 10^ and so the kink does 
not appear. 

In order to study the crossover, we plot in Fig. [5] the 
transmission lengths of H-stacks with N — 10^ and 10'' 
over the wavelength range extended up to A ~ 10^ . As 
for the M-stack, the simulation results for H-stacks can- 
not be distinguished from those of the WSA l^^ . The 
transition from the localized to the near ballistic regime 
occurs without any change in the analytical dependence 
of transmission length, in complete agreement with the 
results of Sec. IIII B 31 The crossover from the near to the 
far ballistic regime is accompanied by a change in the an- 
alytical dependence that occurs at A = X2{N), which for 
these stacks is of the order of lO"' and 10^ respectively. 
The crossover is accompanied by prominent oscillations 
described by Eq. ((5^ . Finally, we note that the vertical 
displacement between the moderately long and extremely 
long wavelength ballistic asymptotes does not depend on 
wavelength, but grows with the size of the stack, accord- 
ing to the law 
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which stems from Eq. ((5(1)) . 

To study the ballistic transmission regime in the re- 
gion A > Ai more closely, we plot in the upper panel 
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FIG. 6; (Color online) Transmission length Zjv vs A for H- 
stacks of = 10^ (solid line), and 10* (dotted line) layers 
(numerical simulation and WSA). Long wave asymptotics for 
the ballistic length in the near and far ballistic regions are 
plotted in thin solid lines. 



of Fig. [Tja) (on a logarithmic scale) the mean value of 
the reflection coefficient for the same two stack sizes. 
That part of each plot which presents the ballistic prop- 
agation comprises two linear asymptotes of the form 
ln(|i?7vP) = const -I- 2 In A, corresponding to different val- 
ues of the constant, one applicable in the near ballistic 
region Ai(A^) ;$ A < A2(A^) and the other in the far 
ballistic region A2(A^) < A. The difference between the 
values of these two constants corresponds precisely with 
the right hand side of Eq. (IM)) , with the far long wave- 
length asymptotic given by the second term in Eq. (|52)) . 
Within the localized regime A < Ai(A^), the reflectance 
is almost unity in all three cases. 

By analogy with M-stacks, the behaviour of the trans- 
mission length in the (A, N) plane is illustrated by the 
phase diagram displayed in the lower panel of Fig. [71 Two 
slanted lines N = /(A) and N = 1{X) divide the plane into 
three parts corresponding to the localization regime (I), 
the crossover regime (II) and the ballistic regime (III). 
The wavelength, at which these lines intersect the hori- 
zontal lines corresponding to the stack lengths, N — 10^, 
and N = 10"*, defines the characteristic wavelengths ((46)) 
for each stack size. From Fig. (TJb), it follows that the 
separation between them grows with increasing size in 
accordance with Eq. (IMl) . It is easy to see that the 
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7: (Color online) (a) Average reflectance from H-stacks 
= 10^ (solid line), and lO** (dotted line) layers (numerical 
simulation and WSA). Long wave asymptotics of the ballistic 
length in the near and far ballistic regions are plotted by thin 
solid lines, (b) Phase diagram for a H-stack. The thick solid 
line corresponds to where the stack size equals the localization 
length. The dashed line corresponds to where the stack size 
equals the crossover length, while the localized, near ballistic 
and far ballistic regimes correspond to regions (I), (II), and 
(III) respectively. 



wavelengths A 1,2 determined from the phase diagram per- 
fectly bound the near ballistic region in both Fig. [6] and 
Fig. Eta). 

We now consider the transmission length for a H-stack 
computed using a single realization. For extremely long 
stacks (N — > 00), the transmission length becomes prac- 
tically non-random, not only within the localized region 
(as is the case for M-stacks) but also in the far ballistic 
region because of the self-averaging nature of the effec- 
tive dielectric constant (|58)) . For less long stacks, how- 
ever, /at (A) also fluctuates in the far ballistic region. To 
demonstrate this, we have plotted in Fig. [8]the transmis- 
sion length Zjv (solid line) and the transmission length 
ZAr(A) for a single randomly chosen realization (dashed 
line). It is evident that, in contrast to the results for 
the M-stack, the H-stack single realization transmission 
length in the near ballistic region is a complicated and ir- 



FIG. 8: (Color online) Transmission lengths Zjv (solid black 
line) and the transmission length for a single realization In 
(dashed blue line) vs A for a H-stack with = 0.25, Qd ~ 0.2 
and A'^ = 10* layers. Each separate point corresponds to a 
particular wavelength with its own realization of a random 
stack. 



regular function, similar to the well known "magneto fin- 
gerprints" of magneto-conductance of a disordered sam- 
ple in the weak localization regime^. In the far ballistic 
region, these fluctuations are moderated since they van- 
ish in the limit as A'^ ^ 00. To support this statement, 
we display in Fig. |5] the set of separate discrete points, 
each of them presenting /at (A) calculated for a different 
randomly chosen realization. 

In summary, we note that the results of Sec. IIV Al show 
excellent agreement between the numerical simulations 
and the analytical predictions of the weak scattering ap- 
proximation of Sec, mil Moreover, even for the case of an 
H-stack of length N — 10'^ and strong disorder, Q^, — 0.9 
and Qd = 0.2, the results of direct simulation and those of 
the WSA analysis coincide completely in the long wave- 
length region and differ by only a few percent in the short 
wave region where the scattering is certainly not weak. 
This is because the perturbation approach based on Eqs. 
(fTTj) . (IT^ is related not to the calculated quantities but 
to the equations they satisfy. 



B. Refractive-index disorder 

1. Suppression of localization 

Here, we present results for stacks with only refractive 
index disorder (RID), with Qi, — 0.25. For H-stacks, 
the transmission length demonstrates qualitatively and 
quantitatively the same behaviour as was observed in 
the presence of both refractive index and thickness dis- 
order. Corresponding formulae for the transmission, lo- 
calization, and ballistic lengths can be obtained from the 
general case by taking the limit as Qd — > 0. 

In the case of M-stacks, however, the situation changes 
markedly. While in the short wavelength region and in 
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the ballistic regime the numerical results are still in a ex- 
cellent agreement with the predictions of the weak scat- 
tering approximation, the WSA fails in the long wave- 
length part of the localization regime. This discrepancy 
manifests itself also for short stacks, in which the localiza- 
tion regime in the long wave region is absent. In Fig.|9]we 
present the transmission length spectrum for a M-stack 
of = lO'^ layers, showing that, for A ~ 5, the numerical 
results for the transmission lengths of M-stacks differ by 
an order of magnitude from those predicted by the WSA 
analysis, and also those observed for the corresponding 
H-stacks. 




A 

FIG. 9: (Color online) Transmission length Zjv vs wavelength 
A for stacks of A = 10^ layers. M-stack: numerical simulation 
(solid line), WSA (dotted middle line). H-stack: numerical 
simulation and WSA (both dashed lines). 

For longer stacks, the differences in the transmission 
length spectra exhibited by M-stacks and H-stacks be- 
comes much more pronounced. In Fig. I10[ we plot trans- 
mission length spectra for different values of the M-stack 
length with A^ = 10^ (using 10^ realizations), 10^ (using 
10'^ realizations) and 10^^ layers (using only a single re- 
alization), with the dashed-dotted straight line in Fig.fTOl 
showing the long wavelength ballistic asymptote ([33| . In 
the moderately long wavelength region corresponding to 
the localization regime, the transmission length In ^ N 
coincides with the localization length I. It exceeds the 
H-stack localization length by a few orders of magnitude 
and is characterized by a completely different wavelength 
dependence. This substantial suppression of localization 
was revealed in our earlier paper— where the localization 
length In was reported to be proportional to A^, in con- 
trast to the classical A^ dependence ([55)) that is observed 
for H-stacks, and which is also valid for M-stacks with 
both refractive index and thickness disorder. The rea- 
son for this difference is the lack of phase accumulation 
caused by the phase cancelation in alternating layers of 
equal thicknesses. 

To study this behaviour in more detail, we generate 
a least squares fit In — AX^ to the transmission length 
data. Respectively, for A^ = 10^, 10^, and lO^'^ layers. 




FIG. 10: (Color online) Transmission length In vs wavelength 
A for an M-stack with = 0.25, Qd = and A = lO'' 
(dashed line), A = lO'' (solid line), and A = 10^^ (dotted 
line) layers. 

the best fits are In = 4A*^■^^ In = OASX"^-^^ and Zjv = 
O.OIA^'^*, indicating that the asymptotic form for the 
localization length differs from a pure power law, and 
perhaps is described by a non-analytic dependence. 

In the crossover part of the long-wave region where 
In ~ A^, the transmission length of M-stacks also differs 
essentially from that for H-stacks. Moreover, the width of 
the crossover region on the Z^v— axis remains of the same 
order of magnitude, although on the A— axis it grows for 
A^ = 10^ to four orders of magnitude, being much wider 
for longer stacks. In the long wavelength region corre- 
sponding to the ballistic regime, the transmission length 
of the RID M-stack coincides with that of RID H-stack. 



2. Transmission resonances 

An important signature of the localization regime is 
the presence of transmission resonances (see, for exam- 
ple, Refsi^"— ), which appears in sufficiently long, open 
systems and which are a "fingerprint" of a given real- 
ization of disorder. These resonances are responsible for 
the difference between two quantities that characterize 
the transmission, namely (ln|Tp) and ln(|rp). The for- 
mer reflects the properties of a typical realization, while 
the main contribution to the latter is generated by a 
small number of almost transparent realizations associ- 
ated with the transmission resonances. 

The natural characteristic of the transmission reso- 
nances is the ratio of the two quantities mentioned above: 

_ (in|rp) 

ln(|TP) • 

In the absence of resonances, this value is close to unity, 
while in the localization regime s > 1. In particular. 
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FIG. 11: (Color online) Ratio s(A) vs. wavelength A for 
— 0.25 and the stack length TV — 10"^. Solid and dashed 
curves are for the RID H-stack and H-stack with Qd ~ 0.2, 
respectively. The middle dashed-dotted curve is for an M- 
stack with Qd ~ 0.25, and the bottom dotted line is for a 
RID M-stack. 



in the high-energy part of the spectrum of a disordered 
system with Gaussian white-noise potential, this ratio 
takes the value 43. 
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FIG. 12: (Color online) Single realization transmittance |Tp 
vs wavelength A for RID M-stacks with = 0.25 and Qd = 



for iV 
line). 



10 layers (solid line) and A'' = 10 layers (dotted 



In Fig. [TTl we plot the ratio s(A) as a function of the 
wavelength for RID M- and H-stacks and for the corre- 
sponding stacks with thickness disorder. In all cases, the 
stack length is iV = 10'^ and it is evident that for the 
RID M-stack s(A) w 1, i.e.. the stack length is too short 
for the localization regime to be realized. In other three 
cases, however, s(A) > 2, which means that the localiza- 
tion takes place even in a comparatively short stack. 

Thus, there are two ways in which to introduce trans- 
mission resonances. The first is to increase the length 
of the stack. Fig. [T^] displays the RID M-stack transmit- 
tance |rp for a single realization as a function of A for two 



FIG. 13: (Color online) Single realization transmittance |rp 
vs A for M-stack of iV = 10^ layers with Q„ = 0.25. Solid line 
corresponds to an M-stack with Qd = 0.2, and the dashed line 
to M-stack with no thickness disorder, i.e., Qd — 0.0. 



lengths: N = 10^ (sohd hue) and N = 10^ (dotted line). 
It is readily seen that while there are no resonances in 
the shorter stacks, they do appear for the longer sample. 
The second way to generate transmission resonances is 
to introduce thickness disorder. To demonstrate this, we 
plot in Fig. [T3] the transmittance of a single M-stack with 
both thickness and refractive index disorder. It clearly 
shows that while the RID M-stack is too short to exhibit 
transmission resonances at A > 3, resonances do emerge 
at longer wavelengths for the M-stack with thickness dis- 
order. 



Effects of the thickness disorder and uncorrelated paring 



10' 
10^ 
10'' 
10^ 
10^ 
10^ 
10^ 



10. 




10° 



10' 
A 



10^ 



FIG. 14: (Color online) Localization length I vs. wavelength 
A for a M-stack with Q^ = 0.25 and Qd=0.001, 0.005, and 
0.01 (from top to bottom). 

Here, we analyze the effect of thickness disorder on the 
A^-anomaly — that is the A^ dependence of the transmis- 



17 



sion length. In Fig. [Ml we plot the transmission length 
In for an M-stack with fixed refractive index disorder 
(Qiy — 0.25) for various values of the thickness disorder. 
It is evident that the transmission length changes from 
Z oc to the classical dependence ^ cx as the thick- 
ness disorder increases from Qd = 0.001 (top curve) to 
Qd = 0.01 (bottom curve). In all cases, the number of 
layers TV = 10^ is longer than the transmission length, 
guaranteeing that Fig. [14] represents the genuine local- 
ization length I. 




10' 



FIG. 15: (Color online) Transmission length In vs. wave- 
length A for the H-stack with R-layers (dashed line) and an 
M-stack (solid line) in which each subsequent layer is chosen 
with equal probability to be of R- or L-type. The stacks in 
both calculations are of the same size. A'' = 10*. 

We have also found that the anomalous dependence 
I (X A^ (or with higher power) is extremely sensitive to the 
alternation of left- and right-handed layers. To demon- 
strate this, we consider an M-stack of length N = 10^, in 
which each subsequent layer is chosen with equal proba- 
bility to be either right- or left-handed. Figure [T5] shows 
the transmission length spectrum for this case, which 
is almost the same for both the H-stacks and M-stacks. 
The only difference is a fairly modest suppression of lo- 
calization, which occurs within the wavelength interval 
0.5 < A < 2.5. This result confirms that it is the addi- 
tional correlation between left-handed and right-handed 
layers in the alternating stack which is responsible for the 
suppression of localization. 



WSA theory in the limits of short or long waves leads to 
the well-known formula 



1 

In 



1 

^abs 



where In is the disorder-induced transmission length in 
the absence of absorption, and the absorption length is 



^abs ^ 



A 
2t^' 



(65) 



For short wavelengths, is a constant given by Eqs 
I, and (|36p . while for long wavelengths, in either 
the localized or ballistic regimes, its contribution is pro- 
portional to and thus is negligibly small in compar- 
ison with the contribution due to losses, which is always 
proportional to A~^. Accordingly, at both short and long 
wavelengths, the attenuation length coincides with the 
absorption length ([65)) . with disorder contributing signif- 
icantly to the attenuation length only in some interme- 
diate wavelength region, provided that the absorption is 
sufficiently small. 




FIG. 16: (Color online) Attenuation length Zatt and absorp- 
tion length iabs vs. wavelength A for an M-stack with disor- 
der = 0.25, Qd = 0.2 and length N = lO". The upper 
solid line displays the (identical) simulation and WSA results 
for a = 10~*; the lower solid line presents numerical results 
while the dashed line displays WSA results for the same ab- 
sorption value. Absorption lengths H65p for both a — 10~* 
and a = 10~^ are shown by dotted straight lines. 



C. Effects of losses 

In this section, we study the transmission through lay- 
ered media with absorption, which is characteristic of 
real metamaterials. In this case, the exponential de- 
cay of the field is due to both Anderson localization and 
absorption- Jj'^-, and in some limiting cases, it is possible 
to distinguish between these contributions. 

For an M-stack with weak fiuctuations of the refractive 
index, weak thickness disorder and weak absorbtion, the 



The results of the numerical calculations shown in 
Fig. [TCj completely confirm the theoretical predictions 
presented above. For weak absorption a = 10^"*, the 
direct simulation and WSA theory give exactly the same 
result (solid curve in Fig. I16p . Over a reasonably wide 
wavelength range, 10^^ ^ 10'^? disorder contributes 
significantly to the attenuation. For such a stack, the 
characteristic wavelengths are Ai(A^) = 47 and A2(A^) = 
265, implying that the contribution of disorder is signifi- 
cant in all regions, from the short wavelength part of the 
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localized regime to the long wavelength ballistic region. 
Due to the losses, however, there are fewer oscillations 
evident in the transmission length spectrum than in the 
lossless case of Fig. O 

For stronger absorption, a ~ 10^^, the wavelength 
range, over which disorder contributes to the attenuation 
is reduced, as well as the relative value of the contribution 
itself. The agreement between the numerical simulations 
and the WSA calculations is reasonable. 



V. CONCLUSIONS 

We have studied the transmission and localization of 
classical waves in one-dimensional disordered structures 
composed of alternating layers of left- and right-handed 
materials (M-stacks) and have compared this to the 
transport in homogeneous structures composed of differ- 
ent layers of the same material (H-stacks). For weakly 
scattering layers and general disorder, where both refrac- 
tive index and thickness of each layer is random, we have 
developed an effective analytical approach, which has en- 
abled us to calculate the transmission length over a wide 
range of input parameters and to describe transmission 
through M- and H-stacks in a unified way. All theoretical 
predictions are in excellent agreement with the results of 
direct numerical simulations. 

There are remarkable distinctions between the trans- 
mission and localization properties of the two types of 
stacks. When both types of disorder (refractive index 
and layer thickness) are present, the transmission length 
of a H-stack in the localized regime coincides with the 
reciprocal of the Lyapunov exponent, while for M-stacks 
these two quantities differ by a numerical prefactor. This 
is quite surprising and, to the best of our knowledge, it is 
the first time that a system where such a difference exists 
has been discovered. 

It is shown that the stacks of M- and H-type mani- 
fest quite different behaviour of the transmission length 
as a function of wavelength. This difference is most 
pronounced in the long wavelength region. In the lo- 
calized regime, stacks of both types are strongly disor- 
dered and reflect the incident wave almost entirely. In 
the ballistic regime, where stacks are almost transpar- 
ent, the transport properties of H-stacks and M-stacks 
are markedly different. H-stacks, over the moderately 
long wavelength ballistic region, and M-stacks, over the 
entire long wavelength region, are weakly scattering dis- 
ordered structures. In the extremely long wavelength 
ballistic region, the H-stack becomes effectively uniform 
— a regime which is absent for M-stacks because of the 
intrinsic non-uniformity caused by the alternating nature 
of the structure. The crossover regions between different 
regimes are comparatively narrow. 

The transmission length for a single realization is non- 
random in the localized regime for both types of stacks. 
It fluctuates strongly over the entire ballistic region for 
M-stacks and in the near ballistic regime for H-stacks. 



In the far long wave region for H-stacks, the fluctuations 
apparent in the transmission length are moderate and 
decrease with increasing stack length. In the case of M- 
stacks, the transition from the localized to the ballistic 
regime is accompanied by a change in the wavelength 
dependence of the transmission length. In contrast, for 
H-stacks, the corresponding change occurs in the vicin- 
ity of the transition from the near to the far ballistic 
regime. Again, the crossover regions between the differ- 
ent regimes are comparatively narrow. 

In M-stacks with only refractive-index disorder, Ander- 
son localization is substantially suppressed, and the lo- 
calization length grows with increasing wavelength much 
faster than the classical square law dependence. The 
crossover region becomes significantly wider, and trans- 
mission resonances occur in much longer stacks than in 
the corresponding H-stacks. 

The effects of absorption on the one-dimensional trans- 
port and localization have also been studied, both an- 
alytically and numerically. In particular, it has been 
shown that the crossover region is particularly sensitive 
to losses, so that even small absorption noticeably sup- 
presses the oscillations of the transmission length in the 
frequency domain. 
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Appendix: Uniform distribution of fluctuations 

If the fluctuations ^j"'' and (5^'''' are uniformly dis- 
tributed in the intervals [—Qi,, Qv] and [—Qd, Qd\ respec- 
tively, and absorption is the same in all slabs (ctj — a), 
the transmission, localization, and ballistic lengths given 
by Eqs. ((22l) - ((26|) can be calculated explicitly in the weak 
scattering approximation. The results are: 



1 



(A.l) 



piA+(l-|-iCT) 

SikQd 



sin(A+Q^) sin(A_Q^)) 



l + ia 



A_g, 



it' 



(A.2) 
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Here, 



+ 2crlm(r) - 2Re(iJi + H2). (A.3) 



where H2 is obtained from Hi by the replacement of A+ 
by A_, and Ei{z) is the exponential integral given by 



[Ei{ilS,+n^) - Ei{iA+n+)] 



00 g-t 



A± = 2k{l±Qd), Et{z)^-J —dt. 

n± — 1 ± + ia, 

SkQ^QdA 

A+Q^cos(A+Q^)] 
i(l + icr^2 
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